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THE ELECTROMAGNETIC TORQUES ON SPHERICAL EARTH SATELLITES 


IN A RAREFIED PARTIALLY IONIZED ATMOSPHERE 


By Frank Hohl 
Langley Research Center 


SUMMARY 


A theoretical investigation has been made of the torques acting, as a 
result of electromagnetic interaction, on spherical conducting earth satel- 
lites. The analysis has been applied to the ^1 -meter -diameter satellite 
Echo II (1964 4A), which is in a near polar circular orbit. The calculations 
have yielded quantitative values for the accelerating and the decelerating 
torques. A probable explanation of the nearly constant spin rate of Echo II 
for a long period of time is that the eddy-current torque and the induction 
torque were in balance. Other torques are shown to be negligible in comparison 
with these two. 


INTRODUCTION 


For a number of scientific experiments aboard artificial earth satellites, 
a knowledge of the satellite motion about its center of mass is required. The 
determination of the satellite motion is a complicated problem and many factors 
must be accounted for. Especially, the torques acting on the satellite must be 
known. It is also of interest to calculate the torques on a satellite to 
determine whether these torques are large enough to explain the fluctuation in 
the spin rate observed for some satellites - for example, the change of the spin 
rate of Sputnik I (1957 Alpha 2) reported in reference 1. The present paper 
reports a theoretical investigation of the torques acting on a spherical arti- 
ficial earth satellite as a result of electromagnetic interaction with the 
ionized atmosphere and the earth's magnetic field. Apparently, some of the 
sources of torque have not been considered heretofore. In the course of the 
calculations, the currents flowing in the satellite hull are also determined. 

A knowledge of these currents may be desirable since they can influence certain 
scientific experiments aboard the satellite. 

The present theory Is applied to the ^1-meter-diameter satellite Echo II 
( 196 ^- 4a) . Calculation of and numerical results for the various torques acting 
on this satellite are presented. One peculiar aspect of Echo II is that the 
period of rotation of this satellite remained very nearly constant at about 
100 seconds for a period of 280 days after launch. However, it is well known 
that the spin rate of the satellite as It spins in the earth's magnetic field 



should he reduced hy eddy-current damping with a time constant of only a few 
months. Thus, one is forced to look for accelerating torques to explain the 
nearly constant spin rate of Echo II. The initial spin-up of the satellite 
cannot he attributed to the torques calculated in this report. It is thought 
to he due to escaping gas from a prohahle rent in the satellite during and 
shortly after inflation. 

The information on the construction and the spin rate of Echo II was 
obtained from an unpublished post launch account and analysis of the perform- 
ance of Echo II prepared hy the Goddard Space Flight Center. 
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The mksAK system of units is used herein, 
area 

magnetic induction 
charge on proton 

unit vector in x-, y-, and z-direction 

electric current 

electric current density 

moment of inertia 

Boltzmann constant 

mass of particle 

magnetic dipole moment 

number density 

power 

charge on satellite 

spherical coordinates 

resistance 

surface resistivity 

time 

torque or temperature 
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velocity- 

Cartesian coordinates 
axes 

angle defined by equation ( 3) 
satellite colatitude 
current flux function 

change in effective satellite radius defined by equation (9) 
angle defined in figure 10 
quantity defined by equation (H) 
free space dielectric constant 
quantity defined by equation (6l) 
quantity defined by equation ( 6 3 ) 
quantity defined by equation ( 7 ) 

radial distance from center of satellite in Debye lengths 
angle between satellite spin axis and magnetic field 

Debye length, 

permeability of free space 

convergence parameter defined by equation ( 50 ) 

mass density 

surface charge density 

inclination of current loop to XY -plane (see fig. 13) 
electric potential 

satellite potential on a circumference where fv s X b) • ? s = 0 
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$ magnetic flux 

cjj angular velocity 

Subscripts : 
a ambient 

e electron 

eff effective 

E earth 

i ion 

m summation index 

n normal component or summation index 

p photoemission 

s satellite 

x,y,z x,y, and z components 

9,i(r 0 and i|r components 

Mathematical notation: 

-> vector 

unit vector 

( ) average over one earth orbit 

| | absolute value 

SATELLITE ENVIRONMENT AND CONSTRUCTION OF ECHO II 

Echo II (196^- hk) is at an average altitude of about 1200 km and is in a 
near circular polar earth orbit. Table I shows some of the pertinent parameters 
of the atmosphere for an altitude of 1200 kilometers. Atmospheric density, ion 
density, and temperature have been taken from Hanson (refs. 2 and 3) Johnson 
(refs, h and 5). It is difficult to assign a definite temperature to the elec- 
trons and ions at the altitude of interest for the present calculations. 
Satellite data indicate that the temperature and density at high altitudes vary 
considerably with the 11 -year solar cycle and also with the time of day. Since 
the res ul ts are to be compared with satellite data obtained near the minimum of 
the 11 -year solar cycle, an ion temperature of 1000° K has been chosen. There 


h 



is evidence that the electron temperature is higher than the neutral gas tem- 
perature (refs. 3 , 6 , and 7 )* For this reason an electron temperature of 
1600° K was used in the calculations. 

Assumptions underlying the present investigation are that the satellite 
consists of a conducting nonferromagnetic spherical shell and that the satellite 
dimensions are large compared with a Debye length. By using the equations given 
"by Spitzer (ref. 8 ) to determine the collision times for ions and electrons, the 
mean free' paths for ions and electrons at an altitude of 1200 kilometers are 
shown in table I to be much larger than the satellite dimensions. Also, the 
number of revolutions about the magnetic field lines between collisions is 
large for both electrons and ions. Thus, the assumptions that there are no 
collisions and that the charged particles move along the magnetic field lines 
are justified. Another assumption is that the ions are singly ionized so that 
n.j^ a = n e ^ a . Experimental results by Bourdeau, Whipple, Donley, and Bauer 

(ref. 9) obtained from the Explorer VIII (i 960 Xi l) satellite support the 
assumption used herein that the electrons and the ions have Maxwellian velocity 
distributions. 

Since a satellite cannot sustain a net current flow to or from the ionized 
atmosphere, the net flux of ions to the satellite must equal the electron flux. 
Since the thermal velocity of the electrons is much greater than that of the 
ions, the satellite will acquire a negative potential. If there is consider- 
able photoemission, however, the satellite may acquire a positive potential. 

The method for calculating the ion and electron accretion by the satellite is 
similar to the method described by Hohl and Wood in reference 10. 

Some of the parameters pertaining to Echo II are presented in table I. 

Echo II is a spherical shell of alum i num foil k.G X 10 "6 meter thick bonded to 
both sides of Mylar polyester film 8.9 X 10 meter thick. There is no elec- 
trical continuity between the parts of the inner aluminum layer. By using the 
thickness of the outside aluminum layer of the satellite and the resistivity of 

aluminum at 313° K (3-1 X 10""6 ohm-meter), the surface resistivity S is found 
to be 6.7 X 10~3 ohm. 

The radius of Echo II is 20.5 meters and the mass of the satellite is 
252.9 kilograms. The satellite rotates at a rate of about O.O 63 radian per 
second. Data received from two beacons located 180° apart on the equator of 
Echo II indicate that the spin axis of the satellite is relatively fixed with 
respect to both space and the satellite. Evaluation of the data at the Goddard 
Space Flight Center also indicates that the spin axis is roughly parallel to 
the earth's magnetic dipole axis. With respect to the satellite, the spin axis 
is roughly normal to a plane which contains the line connecting the two beacons 
at the satellite equator. Thus, it is not determined whether the satellite 
spin axis points through the satellite pole caps, or lies in a plane containing 
the satellite equator, or lies anywhere between these two extremes. 
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QUALITATIVE DESCRIPTION OF ELECTROMAGNETIC TORQUES 


For the purpose of calculating the electromagnetic torques, the earth’s 
magnetic field can he approximated by a dipole field. However, over the 
dimensions of a satellite the earth's magnetic field can be taken to be plane 
parallel in a direction determined by the direction of the earth's dipole field. 
The approximation of the earth's magnetic field by a dipole field can be used 
until effects of the solar wind destroy the dipole field. The effects of the 
solar wind are predominant at altitudes over about 9 earth radii. Thus, the 
dipole approximation should be valid out to several earth radii. 

The electromagnetic torques acting on a satellite can be divided into two 
types. Torques of the first type result from the interaction of the metallic 
satellite hull with the earth's magnetic field. The relative motion of the 
satellite with respect to the earth's magnetic field may induce currents in the 
satellite hull; these currents can then interact with the earth's magnetic 
field to produce torques. Also, as the satellite passes across lines of force 
of the earth's magnetic field, a polarization of the satellite occurs. This 
polarization results in a nonuniform surface charge. If the satellite now 
rotates about its spin axis, the surface charge will keep its orientation with 
respect to the magnetic field and the satellite velocity vector and thus con- 
stitute a motion of electric charge relative to the conducting shell. The 
Joule heat dissipated by the current can be shown to come from the energy of 
rotation of the satellite. Thus, the first type of torque represents a damping 
torque acting on the satellite and would be present even if the satellite were 
moving in a nonionized atmosphere. 

Torques of the second type result from the interaction of the satellite 
with the ionized atmosphere in the presence of the earth's magnetic field. 

These torques are the Coulomb torque and the induction torque. The Coulomb 
torque is due to unsymmetrical accretion of electrons and ions. If the satel- 
lite velocity is parallel to the earth's magnetic field or if there is no mag- 
netic field, the satellite is unpolarized and the satellite cross section for 
ion impingement is as shown in figure 1. As indicated, the negative charge on 
the satellite causes the cross section for ion impingement to be increased 
symmetrically so that no Coulomb torques are developed. If, however, the satel- 
lite moves across lines of force of the earth's magnetic field, it becomes 
polarized and the cross section for ion impingement is increased on the more 
negative side of the satellite and decreased on the more positive side, as 
shown in figure 2. The unsymmetrical impingement of ions produces a torque. 

As was indicated by Hohl and Wood in reference 10, the momentum transferred by 
the ions that are scattered but that do not impact the satellite is negligible 
compared with the momentum transferred by the ions that impact the satellite 
as a result of the increase in effective cross section due to the charge on the 
satellite. It should also be noted that electrons contribute a negligible 
amount to the Coulomb torque, since the electrons impinge on the satellite from 
two sides along the magnetic field lines and the electron mass is much smaller 
than the ion mass. 
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The induction torque is the result of currents flowing in the satellite 
hull due to preferred accretion of electrons or ions on certain portions of the 
satellite surface. The induction torque is illustrated in figure 3* Ions 
impinge primarily on the front portion of the satellite, "being concentrated 
somewhat on the more negative side of the satellite because of the increased 
cross section for ion interception. Electrons moving 'along the magnetic field 
lines impinge on the satellite principally on the side which is more positive. 

Qualitatively, the current is then as indicated in figure 3j in general, the 
— > — ) 

resulting j X B force does not pass through the center of the satellite and 
therefore it exerts a torque on the satellite. The energy put into rotation is 
taken from the kinetic energy associated with the satellite motion around the 
earth. 

A nonelectrodynamic torque of interest is the aerodynamic torque due to 
the satellite spin. It is estimated herein in order to he compared with the 
other torques. In the calculation done in this report the satellite spin vec- 
tor is pointing in the negative z-direction shown in figures 4 and 5* 


SATELLITE POTENTIAL 


In the absence of photoemission, an artificial earth satellite will 
acquire a net negative charge. The reason for the acquisition of the negative 
charge is that for an uncharged body moving at satellite velocity the rate of 
accretion of electrons is much greater than that of ions because the thermal 
velocity of the electrons is much higher than the velocity of the satellite or 
the thermal velocity of the ions. To keep the rate of electron accretion equal 
to that of ion accretion, the satellite will acquire a negative potential. 

The equilibrium potential of the satellite is obtained by setting the sum 
of the electron current and the ion current to the satellite equal to zero - 
that is. 


ii 


+ i e - 0 


( 1 ) 


It is assumed that upon striking the satellite surface the ions become neutral- 
ized and either stick to the surface or are reflected or reemitted as neutral 
particles. 

Experimental results of Hagstrum (refs. 11, 12, and 13) show that ions 
striking a metal surface will be neutralized by Auger charge exchange with an 
efficiency of nearly 100 percent. The experiments by Hagstrum also indicate 
that when the ions strike the metal surface some secondary electrons are pro- 
duced. In the calculations that follow, these secondary electrons are 
neglected. 


To determine the ion and the electron currents to the satellite, it is 
necessary to specify the orientations of the magnetic field, the satellite 
velocity vector, and the coordinate system to be used. This information is 
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contained in figure The earth coordinate system is fixed with its origin at 
the earth's center. The satellite coordinate system has its origin at the 
center of the spherical satellite and the Z s and X s axes remain always 
parallel to the Z-g and Xg axes, respectively. Since the satellite is t ak en 
to he in a circular polar earth orbit at an altitude of 1200 kilometers, the 
earth's dipole field in the satellite coordinate system, as determined from the 
information in reference l4, is 


B = B 


a 


-i3 cos p sin p + k ^1 - 3 cos^p^ 


( 2 ) 


with B a = 1.9 X 10"5 weber per square meter. The angle a shown in figure 4 
is given by 


3 cos p sin p 
1 - 3 cos2p 


(3) 


The conditions for Echo II are such that the inequalities V e » V s and 
V s » Vj_ can be assumed to hold. In reference 15, the second inequality is 
shown to result in an error in ion acquisition rate of only 4 percent on the 
front and 4 percent on the rear of the satellite at an altitude of 1200 kilo- 
meters. The kinetic-theory differential expression for the flux of particles 
along the magnetic field lines in one direction (with velocities along the 
field lines between V and V + dV) is 


dj e 


-en^ 


me 


2jtkT ) 


\l/2 


V exp 



dV 


( 4 ) 


To obtain the electron current density to the satellite, equation (4) must be 
integrated between the kinetic energy limits -e0 g for < 0 and °° and 
between the limits 0 and <» for 0 S ^ 0. 

The electrons reach the satellite mainly along the magnetic field lines. 
When referred to the satellite orientation shown in figure 5, 'the electron cur- 
rent density to the spherical surface is 



where the trigonometric factor comes from the dot product of the unit vec- 

/v A 

tors r and B and where 


0s = 0eq + (Vs X B) • ? B 
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0 S = 0 e ^ + 2V s B a r s cos p sin 0 sin i|r ( 6 ) 

and 

e 0 = 1 ( 0 s 

€ 0 = 0 (0 S 

In equation ( 6 ) the ohmic potential drop has been neglected in comparison with 
the induced potential gradient. That the neglect of the ohmic drop is valid is 
seen in a subsequent section when the currents flowing in the satellite hull 
have been determined. 

The symbol €0 is introduced to limit the maximum electron current that 
can be intercepted by the surface to that for which 0 S = 0. Increasing the 
potential 0 S above zero will not appreciably increase the electron current. 

The general behavior of the electron current as given by equation ( 5 ) has been 
experimentally verified by Bourdeau, Donley, Serbu, and Whipple (ref. l6) 
except that they found a small reduction in electron current on the back side 
of the satellite which is probably due to a negative charge in the wake. 

The ion current to the satellite is obtained from the rate at which the 
ions are swept up by the effective satellite cross section for ion interception 
and is given by 



^i n a e ^s^-eff 


(8) 


Ion current measurements on the Explorer VIII satellite at an altitude of 
1000 kilometers are in good agreement with the assumption that the ions are 
incident only on the front surface of the satellite irrespective of the orien- 
tation of the earth's magnetic field. 


The effective cross-sectional area for ion interception & e ff and the 
equilibrium potential 0 e ^ can be determined by successive approximations as 
follows: With a tentative value for 0 e qy ‘the potential distribution around 

the satellite is determined. Then, the ion trajectories are calculated to give 
A e ff. These values of 0 e ^ and A e ff are then used to determine whether the 

net current to the satellite is zero. However, the result given by Jastrow and 
Pearse (ref. 17 ) for the sheath thickness approximates the results that are 
obtained when the actual trajectories are calculated; that is, any ion which 
approaches the satellite with an impact parameter greater than 


r eff 


r a - 


I 0£ 


2e 


Jl 

kT 


7\ d = 


+ S 


(9) 
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will not be intercepted by the satellite, whereas any ion which approaches the 
satellite with an impact parameter less than r s + 6 will be intercepted. The 
effective radius of the satellite for ion interception is then approximated by 
equation (9)- 

Comparisons of the effective satellite radius as obtained from equation (9) 
and as obtained by computing ion trajectories are made in the appendix. 

When referred to the coordinate system shown in figure the ion current 
density to the satellite becomes 

jq = e n a V s ^ ~ff^ (cos £ sin 9 cos \|r - sin p cos 9)eq (lO) 


where 


€ i = 1 


e i “ 0 


cos p sin 9 cos i|r 

- sin p cos 9 ^ 0, 

cos p sin 9 cos \|f 

- sin p cos 0 < 0 


> (H) 




The equilibrium potential is determined from 



j e dA + 


jq dA = 0 


( 12 ) 


where the integration is over the spherical satellite surface. 

Figure 6 gives the equilibrium potential of a spherical satellite as a 

function of its radius (conditions are the same as those for Echo II) for a 

— ^ ^ __ 

typical value of V s X B. The equilibrium potential Is a function of V s , B, 

and r g and decreases nearly linearly with increasing r s . If equation (12) 

is written explicitly, it can be deduced that the equilibrium potential will 
decrease nearly linearly with increasing B. 

In figure 7(a) the equilibrium potential is given as a function of colati- 
tude for Echo II. The results of figure 7(a) Indicate that for Echo II the 
equilibrium potential can be approximated by 


0 e q = -0.2 - 4.56 |cos pj 


(15) 
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and, therefore, the satellite potential can he approximated by 


0 S = -0.2 - 4.561003 p| + 5*67 cos p sin 0 (l4) 

* 

Experimental results show that the effect of photoelectric emission is 
generally not negligible; therefore, photoemission should he included in the 
calculations. Chang and Smith (ref. 18) and Brundin (ref. 19 ) have calculated 
the photoelectric emission from a satellite hy considering the radiation inci- 
dent on the satellite and hy using an appropriate absorption coefficient and 
work function. However, since experimental measurements of photoemission are 
available, the measured currents can he used to determine the effects of photo- 
emission on the equilibrium potential and on the torques. Considering the 
results obtained by Bourdeau, Donley, Serbu, and Whipple (ref. l6) from 
Explorer VIII data, a reasonable value for the photoelectric current is 
5 X 10 “5 ampere per square meter. If the satellite is in sunlight, the equi- 
librium potential is obtained from 




dA + 5 X 10"5 rt r s 2 = 0 


(15) 


The resulting equilibrium potential as a function of colatitude is given in 
figure 7("b)* The equilibrium potential with photoemission can then be written 
as 


0 e q = - 0.165 - 1.688 |cos p| (l6) 

If the satellite surface has a positive potential, the photoelectric emission 
current will be reduced. The dependence of the photoelectric emission current 
on the positive satellite potential has been neglected in the foregoing calcu- 
lations. Since most of the satellite surface is always negative, the error 
introduced should not be large. 


RETARDING TORQUE DUE TO SURFACE CHARGE 

— > — > 

The induced potential gradient on the satellite is V s X B. Thus, there 
exists a potential difference across the satellite which equals the induced 
potential difference minus the ohmic voltage drop; this ohmic voltage drop, 
however, is negligible. Therefore, the potential difference across the satel- 
lite is given by 


0 = ^ r s v s B a cos P 


(17) 
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The maximum potential induced across Echo II (over the earth's poles) is 
found to "be about 11 volts. The potential difference polarizes the satellite 
by producing a surface charge. As the satellite rotates, the orientation of 
the surface charge with respect to the earth's magnetic field and the satellite 
velocity will remain the same Thus, there is an electric current over the 
satellite surface. The resulting Joule losses tend to reduce the 'rate of rota- 
tion of the satellite. 

To obtain the surface charge, the potential distribution through the 
sheath could be calculated numerically. However, if it is assumed that the 
satellite is in vacuum, the variation of potential around the satellite is 
given by a solution (ref. 20) of Laplace's equation - namely, 

0eq r s __ V s Br s 2 
r “ r 2 

where the orientation of the satellite coordinate system has been so chosen 
— ^ — ) 

that V s is along the Y-axis and B is along the X-axis. Inasmuch as the 
surface charge density a is equal to the normal component of the dielectric 
displacement at the surface. 



a = 



h = e Q ” - 2e 0 V s B cos 0 (19) 

s 

To determine the maximum torque, the satellite is assumed to rotate about 
the Y-axis, because the resulting flow rate of surface charge is then a maximum. 
In the coordinate system in which the surface charge is at rest, the surface 
current density due to the satellite rotation is 


— » 

j 


(r s Xffl) 


j = 2e 0 V s Br s U) cos ©( i cos 0 - k sin 0 cos \|r. 


( 20 ) 


The rate at which the current liberates Joule heat is 


P = 


j 2 s dA = ^n(e o V s r s 2 B03) 2 S 


15 


( 21 ) 


The energy source for the Joule heat is the kinetic energy of the rotating 
shell. Differentiating the equation of kinetic energy for the rotating 
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spherical shell gives 



and thus the corresponding decelerating torque is 


( 22 ) 


t = = -£ 

dt no 


T = -SijtSm 
15 


(e 0 V s r s 2 B f = -6 


X 


10“ 20 cn 


newton-meter 


(23) 


In equation ( 23 ), the parameters pertaining to Echo II have been used. For any 
reasonable value of «), the decelerating torque is completely negligible in 
comparison with other torques that are considered herein. 


RETARDING- TORQUES DUE TO EDDY CURRENTS 


It is well known that in the presence of a magnetic field the angular 
velocity of a conductor tends to be decreased (ref. 21). If the satellite 
rotates about an axis parallel to the magnetic field, no eddy currents are 
induced because of the rotation. If, however, the satellite rotates about an 
axis which has a component perpendicular to the magnetic field, eddy currents 
are induced in the conducting hull of the satellite. A number of papers have 
been written on the subject of eddy-current torques (for example, refs. 22, 23 > 
and 24) . Thus for Echo II only the torque for two specific orientations are 
considered. 


The retarding torque due to the eddy currents in a uniformly conducting 
spherical shell has been shown by Halverson and Cohen (ref. 25) to be, when the 
spin axis is parallel to the Z-axis, 


T = 


5izur s 


V 


9S 2 + p n ^n 2 r 2 


( A A \ A 

-I3S + jp^urg) sin 2r) + k6S sin 2 t 


(24) 


For Echo II S » p. 0 tDr s so that, if the satellite were uniformly conducting, 
the torque would be 


-» jtrnrc i^ 2 / ~ ^ o \ 

T ^ ^-i sin 2tj + k2 sin 2 qj 


(25) 


where q is the angle between the spin axis of the satellite and the magnetic 
field. Since the satellite spin axis is along the Z-axis for the present 
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problem, the angle t) is equal to the angle a indicated in figures 4 and 5* 
The satellite coordinate system shown in figures 4 and 5 is that used for most 
of the calculations. On substitution of equation (2) for B in equation (25) 
the torque becomes 

T = - — (l + 5 cos^p) ^-i sin 2a + k2 sin^a) (26) 

The components of the eddy-current torque for the parameters of Echo II are 
shown in figure 8 as a function of (3. Because of symmetry the torque is shown 
for only one-half orbit. 

This eddy-current torque was calculated -under the assumption that the 
satellite consists of a uniformly conducting spherical shell. However, Echo II 
is made up of 106 individual gores which are electrically connected only near 
the poles, as shown in figure 9- Thus, a better approximation to the torque 
can be obtained by assuming that the current paths are only along the gores. 

To calculate the eddy-current torque for Echo II, the orientation of the 
satellite with respect to the spin axis must also be specified. Two cases of 
satellite orientation with respect to the spin axis are considered: case I - 

the satellite is spinning about an axis through its poles; case II - the satel- 
lite is spinning about an axis which lies in the plane of the satellite equator. 
The satellite orientation for case I is shown in figure 10. Note that the 
orientation of the coordinate system with respect to the earth remains as indi- 
cated in figures 4 and 5 and that B remains in the ZX-plane as indicated in 
figure 10. Consider the eddy currents flowing in circles (along the gores) 
through the poles. If it is assumed that the magnetic field due to the eddy 
currents can be neglected, the magnetic flux through a ring made up of two 
gores on opposite sides of the satellite becomes 

$ = TCr s % sin a sin i|r 


$ = jtr s 2 B sin a s:i - n fJC * 


(27) 


The electromotive force producing current in the ring is thus 


0 = - — = -uurr s ^B sin a cos cut 
dt 


(28) 


The resistance of the ring is 


'rt-e 


E = 2S 


d9 

At| r sin 9 
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where e is the angle indicated in figure 10 to the point where the gores are 
electrically connected to each other and A\|r is the angle corresponding to the 
angular width of the gores. 

The angle e where the gores are connected to the connecting ring is 
0.03 radian. Using in equation ( 29 ) this value of e and 6.7 X 10”5 ohm for 
the surface resistivity S gives the value of the resistance R as 1.8 ohms. 
In this calculation, however, the resistance of the connecting ring and that of 
the connections where the connecting ring is ultrasonically welded to the gores 
were neglected. Measurements on the backup satellite for Echo II indicate 
that the average value of the resistance at each connection (not including the 
connections that were open circuits) was about O .65 ohm. If this value is 
assumed to be applicable to the satellite that was launched, then the value of 
R to be used for the present calculations should be increased to 4-. 4 ohms. 

The current flowing in a ring consisting of two gores is 

1q = & = -§-r s ^Bco sin a cos ^ ( 30 ) 

Ja I\ 

The magnetic moment of the current is 

— ) — ) 

M = iA 


The torque on the ring can now be calculated from 


T = M X B 


= sin a cos ^i sin ^ + j cos \|r) X ^i sin a + k cos a) 


T 


= -g!rt2 rs 4g2 s ^ n a CQS cos a cog ^ _ j cos a s q n t|j _ k sin a cos 7 ) 


(51) 


Averaging over \|r = cot and using B 2 = B a 2 ^1 + 3 cos 2 (3^ gives 


T = 

4r 


'^r s ^B a ^ (l + 3 cos 2 p)(-i sin 2a + k2 sin 2 a) ( 32 ) 
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Since there are 53 such rings, the total torque on the satellite is 

T = ” (l + 3 cos2p)(-i sin 2a + k2 sin^a) (53) 

The torque given by equation (33) would be the correct torque if all the 
gores were electrically connected to the connecting ring. However, electrical 
continuity measurements on the backup satellite for Echo II showed that about 
one -fourth of the connections between the gores and the connecting ring at each 
connecting ring represented open circuits (no continuity) . One can therefore 
assume that about one-half of the gores were not connected either to the top or 
to the bottom connecting ring. The open circuits will very likely decrease the 
eddy-current torque by a factor of more than 2, since both ends of a gore must 
be connected for eddy current to flow. To obtain a better approximation for 
the eddy-current torque experienced by Echo II, the torque given by equa- 
tion (33) is divided by 3* The factor of 3 was arrived at by the author after 
considering probable combinations of open circuits. The resulting components 
of the torque as a function of p are shown in figure 11. It can be seen that 
the torque curves for Echo II are the same as the curves for the uniformly con- 
ducting spherical shell except for a reduction in magnitude. 

The orientation for case II is shown in figure 12. Consider the ring 
inclined at the angle r to the XT -plane. The flux through this ring is 

$ = itr s %(cos a cos r - sin a sin r cos ait) (3^) 

where 

flit = i|c 

Thus, the potential induced in the ring is 

0 = = -niBjrrg^sin a sin r sin ait (35) 

The current flowing along the gore (at t) is 

i = -^n:r s % sin a sin t sin ait ( 36 ) 


Again, the torque on the ring is obtained from 
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— > — > — > 
T = M X B 


= - s ^sin a, sin r sin i|r ^-i sin t cos \(r 


- j sin r sin ^ + k cos r) x (i sin a + k cos a) 


s ^sin a sin t sin [ji sin r cos a sin \|r 
+ j(cos r sin a 4- sin r cos a cos ijr) + k sin r sin a sin (57) 

Averaging over \|r = oat and. using B 2 = B a 2 ^l + 5 cos 2 p) gives 

? =■ -§B a 2 »t2r^(l + } cos%)( -i sin 2 r cos a sin a + k sin 2 r sin 2 a) ( 38 ) 
To get the total torque the sum over t = n is required - that is, 

53 

T = -j^r s \2(l + 3 cos 2 p) (-1 sin 2a + k2 sin 2 a) < y sin 2 (n At) (39) 


The resulting components of this torque as a function of |3 are indicated in 
figure 13 , where equation (39) has been divided by the factor 3 to account for 
open circuits. 

It should be noted that in its orbit around the earth the satellite will 
encounter a varying magnetic field. However, since oi is much larger than the 
rate of revolution of the satellite around the earth, it can be assumed that 
the direction of the field does not change appreciably during a rotation of the 
satellite about its axis. The effect of the time varying earth's magnetic 
field with respect to the coordinate system fixed in the satellite can there- 
fore be neglected. The problem of a satellite which remains stationary rela- 
tive to the coordinate system moving with it in its orbit around the earth has 
been considered by Zonov (ref. 24) . 


COULOMB TORQUE 


The Coulomb torque is due to unsymmetrical accretion of ions. The polari- 
zation of the satellite causes more ions to be intercepted by the more negative 
side of the satellite than by the more positive side. Thus an unsymmetrical 
transfer of momentum occurs which results in a torque on the satellite. The 
contribution to the torque from the ions that are deflected but do not strike 
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the satellite is negligible as indicated by Jastrow and Pearse (ref. 17) and by 
Hohl and Wood (ref. 10) . To obtain the Coulomb torque the ion trajectories 
might first be calculated; this was done by Hohl and Wood in reference 10 for 
Explorer IX (1961 Delta l) . The maximum impact parameters for impingement 
obtained from these trajectories are in reasonably good agreement with the 
theory given by Jastrow and Pearse. Thus, even though the satellite potential 
is not spherically symmetric, the Jastrow-Pearse theory applied to a point of 
the satellite surface gives the approximate thickness of the sheath at that 
point. A comparison of sheath thickness as determined from calculations of ion 
trajectories with that given by the Jastrow and Pearse theory is given in the 
appendix . 

If the Coulomb torque were found to be a considerable fraction of the 
torque required to balance the eddy-current torque, a more exact method would 
have to be used in calculating the Coulomb torque. However, since the Coulomb 
torque will be found negligible compared with some of the other torques, a 
10- to 30-percent uncertainty in it will not affect any of the results. Thus, 
the effective satellite radius for ion interception can be sufficiently well 
approximated by the Jastrow and Pearse result (eq. ( 9 )). 

The calculation of the Coulomb torque is simplified if the satellite 
coordinate system shown in figure l4 is used. The magnetic field B and the 
satellite velocity V s are then in the XY -plane and the satellite potential 0 S 
is given by 

0 s = 0 eq " 2 V s B a r s cos P cos 0 ' (^°) 


If it is assumed that the ions transfer their momentum to the satellite at 
the point where they are intercepted by the sheath, the magnitude of the torque 
can be approximated as 



Equation (4l) gives the torque for the condition in which the ions that strike 
the satellite "stick" to the satellite surface. If the ions strike the satel- 
lite surface and are reemitted, the torque will be different. If the ions 
strike the surface and are reemitted as neutral particles, the contribution to 
the torque resulting from the neutralized ions recoiling from the front surface 
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of the satellite -would have to he considered; that is, most of the ions which 
strike the satellite surface have been accelerated by the electric field of 
the satellite and, hence, if the ions are reflected as neutral particles, 
their final velocity after reflection will be much larger than their initial 
velocity. The radial component of the electric field of the satellite will 
accelerate the ions in such a manner that they cannot contribute to the torque. 
However, since the potential around the satellite is not completely radial, the 
neutralized ions recoiling from the front surface of the satellite with a 
velocity larger than their initial velocity would produce a torque. The compu- 
tation of this torque is complicated because the ion trajectories must all be 
computed. An approximate method used to calculate the torque due to the 
recoiling neutralized ions yields results of the same order of magnitude as 
those obtained by use of equation (4l) for the Coulomb torque. For the present 
calculation, only the torque due to the impingement, and not the reemission, of 
the ions is considered. 

With reference to figure 4, the torque is found to be either parallel or 
anti-parallel to the radius vector from the center of the earth through the 
center of the satellite - that is, 

T(p) = -T(p)^i sin p + k [cos p|^ (42) 

The components of the torque as calculated for Echo II by use of equations (4l) 
and (42) are shown as a function of p in figure 15 • Note that the torques 
are always referred to the coordinate system shown in figures 4 and 5- 


INDUCTION TORQUE 


The induction torque is by far the largest of the accelerating torques and 
is considerably more difficult to compute than any of the other torques. To 
determine the torque resulting from the unsymmetrical accretion of charge, the 
current distribution over the satellite surface must be found. The current 
distribution is calculated for a uniformly conducting spherical satellite and 
for Echo II where the spherical shell is made up of gores that are electrically 
connected together only near the poles. 

To obtain the current distribution for a uniformly conducting spherical 
satellite, the coordinate system shown in figures 4 and 5 is used. By addition 
of equations (5) and (10) the current density to a point on the satellite sur- 
face without photoelectric emission is 

ji + j e = en a^s(~§^^ (cos p sin 0 cos \|r - sin p cos 0)ep 


“ en a 



sin 9 cos sin a + cos a cos 



(45) 
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If the increased satellite cross section for ion interception is neglected in 
the calculations of the induction torque, the resulting error is less than 
1 percent. Using the approximate expression for the effective cross-sectional 
area given by equation ( 9 ) makes the error even smaller. 

— ) 

Let j g he the surface current density over the satellite surface. Then 

— ) 

j s can he expressed as 



m 


where 7 is determined from 

= ~(ji + J e ) (^5) 


In finite difference form equation (45) becomes 




V 2 7 = 


P P 
r s sin 9 


sin 9^- (sin 9^- 

59 V 59 


5 2 7 

5^i|r 


where A9 and At represent the separation of the mesh point used in the 
numerical solution. The values of subscript m are 1, 2,. . ., M and those 
of subscript n are 1, 2, 3,. . . , N. The mesh points are specified by 


fm = ( m 



The boundary conditions are 


- i)At (At = sf) (Vr) 

- |)A9 (A9 = g) (48) 


^W+l,m ^n^M+l ?n,l 


A 


7 °; m 7 l,m± — 




^n,0 ^n,M 


J 


(^9) 
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The equation actually used for the numerical solution on the IBM 709^ data 
processing system was of the form 


■y(r+l) = + v|-pr 

n,m 'n,in 1 ■- 


-1 


V 

L(A6) 2 

sin 9 n cos gn] 

2 A9 

J 

Jr) (1 

J 

7 n,m+l(^ 

*)J 

^m) = "(^e 

+ jq) 


sln29 n + / 
(AQ) 2 V 


-V 

a *) 


<[r s 2 sin 2 0 n f (0 n ^ m ) - 7^^ 


sin 2 0 


n 


(A0) : 


_ J r+1 Vi_\ 2 _ 7 ( r ) 

/ n,m-lyty ) n+l.,m 


sin 2 0 n sin 0n cos^ 9 n 


(A0) : 


2 A0 


7 


n,m 


(50) 


to speed convergence and the magnitude of which was determined hy the rate of 
convergence. The superscript (r) indicates the rth iteration. The torque 
due to the surface currents is determined from 


T = 


r- 2 it n rt 

0 0 


[r s x(? B X B ) r s 2 sin 0 


d0 d\|r 


or 


T x = r 8 5B a (l + 3 cos 2 p) l/2 


j' 2lt n rt 

0 ^ 0 

+ cos 0 cos cos a - sin 9 sin 2 \(r sin a) 


|Jq cos 0(sin 0 sin a 


(51) 


- sin >|r(sin 9 cos \|r sin a + cos 9 cos a)J sin 9 d0 d\[r ( 52 ) 


ly - r s' 


/ \-l / 

p 2it p ^ t 

( 1 + 5 cos 2 pj ' 

\ / [ 
J Q W Q 


+ sin 0 cos \|r sin a) + cos i|r(cos 9 cos a 
+ sin 0 cos i|r sin a)j sin 9 d9 dip 


(53) 


21 



%a ^1 + 3 cos 2 p J" J' [^0 s ^- n ®(sin 9 cos y sin a 


T z = r s 


- cos 0 cos a^j sin 0 &0 d\|r 

The two components of j* are given by 

^n-l.m “ ^n+l.m /. \ 
j0 " 2 A0r s = Wn, 


(5^) 


(55) 


and 


- ^n^-i- ~ ^n,m+l /. \ 

^ 2 Afr s sin 0 ' ^'n,m 


(56) 


By using the results of the numerical solutions for Jq and the compo- 

nents of T are obtained from 


T x = r s 


5B a (l + 3 cos 2 p ) l/2 XI («Je)n m cos 9 n ( sln 0 n sin 

n m x 

+ cos 9 n cos t m cos a - sin 0 n sin^i)f n sin 
- sin (sin 0 n cos t m sin a + cos 0 n cos 


a 


sin 0 n A0 A)r (57) 


T y = r s 


% a (l + 3 COS 2 P ) 1/2 £ ? |je)n ;m cos 9 n sin ^m( cos 9 n cos a 


T = r 
z x s 


+ sin 0 n cos \|r n sin a) + m ( cos 7 m cos 0 n cos a 

+ sin 0 n cos if m sin ajj sin 0 n A0 A»|r ( 58 ) 

[(jol m sin 9 n( sln 9 n cos sin a 
n m t- ' 


- cos 0 n cos 


a) 


sin 0 n A0 A'Sf 


(59) 


The induction torque for the uniformly conducting satellite without photo- 
emission as given by equations (57) and (59) is shown in figure 16. Wo 
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Ty component is shown in figure 1 6 , since the computed value of Ty was 
smaller than the error in the numerical computations. 

If photoemission is included in the determination of the induction torque, 
then the orientation of the satellite orbit with respect to the earth-sun line 
must be specified. The induction torque with photoemission is calculated for 
two orientations. With reference to figures 4 and 5, for orientation I the 
solar radiation is taken to be incident along the X-axis of the satellite 
coordinate system. Thus the torque needs to be calculated only for the part of 
the orbit during which the satellite is in sunlight. The function f (9n>^m) 
appearing in equation ( 50 ) must be modified from the expression given in equa- 
tion (43) ty adding the photoelectric emission current, as used in equation ( 15 ) 


where 


= 5 X 10 ”5 cos ty sin Se-r, 
P & 


(6o) 





(cos ty ^ 0) 


(cos ty < 



(61) 


The calculation of the current distribution and of the induction torque 
for the satellite with photoemission then proceeds the same as that for the 
uniformly conducting satellite without photoemission. The resulting torque for 
the uniformly conducting satellite with photoemission for orientation I is 
given in figure 17(a). Note that for one-half of the orbit the torque without 
photoemission is shown since during about one-half of the orbit the satellite 
is in the earth's shadow. 

With reference to figures 4 and 5, the satellite orbit for orientation II 
is such that the solar radiation is incident along the Y-axis of the satellite 
coordinate system. The satellite will then remain in sunlight during the 
entire earth orbit. Now, f^n^m) must "be modified from the expression given 
in equation (43) by adding 


= 5 X 10 ”5 S in ty sin 06^ 


( 62 ) 


where 


t 



1 


e 


P 


= 0 


(sin ty 5? 0) 
(sin ty < 0) 


( 63 ) 


The resulting torque for the uniformly conducting satellite with photoemission 
for orientation II is given in figure 17(b). 
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The torque for a spherical satellite made up of 10 6 gores is calculated 
next. It is assumed that the current flow is only along the gores. Also, the 
gores are connected near the poles so that the current can flow from one gore 
into another. For calculation of the torque the orientation of the satellite 
with respect to its spin axis is shown in figure 18. 


The current in each gore must he 
element of surface of the nth gore is 


determined. The current flowing into an 
given hy (see eqs. (5) and (10)) 


dl n (|3,e) 



P sin 9 cos f n 


- sin p cos 0 


h 


- Isin 0 cos 

]j2 Jang I 

+ cos a cos 0I exp 


t n sin a 



Ai|/r s ^sin 0 d0 


(64) 


where 



(n =0,1,2, . . .,105) (65) 


and a, e^, e^, and 0 S are given hy equations (3), (ll), (j) , and- (6), 

respectively. The total current flowing in the nth gore (in the 0-direction) 
at 0 is then 


i n (p,0) = l n (p,0) + I n (p) 

(66) 

In(P>9) = f cLEn(fM) 

( 67 ) 

In(P) = rE )(P) “ 

(68) 

R = [*" 3 cL9 

J e At sin 0 

(69) 

nit-e 

0 n (P) = / l n (p,0)— 4— d0 

J g Ai|r sin 0 

(70) 
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(7D 


105 

^ = isr 71 ^ {|3> 


n=0 


The resistance R as given in equation ( 69 ) is the resistance of one gore 
between the two ends where it is electrically connected to the other gores and 
e is the angle from the polar axis to the connecting ring as shown in fig- 
ure 10. An additional current I n (p) must flow in the nth gores of such mag- 
nitude that the voltage difference across that gore is equal to the average of 
the voltage differences across all gores due to the currents i n (p,9). 

The force on an element of a gore d l (fig. l8) is then given by 

® n =i n (p,0)dUB (72) 

The corresponding torque is 

cLT n = r X dP n (73) 


The resulting torque on one of the gores is 

-> it-e 


T n (p) = r s 2 B a (l + 3 cosSp) 1 / 2 ^ 

+ cos 0 cos <x^i n (9,p)d0 

+ cos a cos 0^i n ((3,9)d0 
+ cos 0 cos a^i n (p,0) 


^n( sl 


cos 9 cos \lr n ( sin 0 cos 7 n sin a 


+ J 


- k 


' it-e 


' jt-e 


cos 9 sin tyn ^sin 9 cos \|/ n sin a 


sin 9 (sin 9 cos i]/ n sin a 


d0 


( 7*0 


By summing over all gores, the total torque becomes 

105 

T(P) = 21 Tn(p) 

n=0 


(75) 


Equation (75) has been numerically integrated. The resulting torque as a func- 
tion of p for Echo II without photoemission is shown in figure 19- 
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As was done for the -uniformly conducting satellite, photoemission is next 
included and the induction torque is calculated for two orientations of the 
satellite orbit with respect to the earth-sun line. Again, for orientation I, 
the solar radiation is taken to be incident along the X-axis of the satellite 
coordinate system. Equation (6^) must then be modified to 


di n (p,e) = 


n a eV s (cos (3 sin 0 cos - sin p cos 
- n a e i s ^ n ® cos ^n s ^- n a + cos a cos ®l 


+ 5 x 10 ~ 5 s in 0 cos \|r n e 


■P 


Ai|rr s ^sin 9 d0 


( 76 ) 


With expression (j6) being used for dl n (p,9), the method of calculating the 
induction torque for Echo II with photoemission is exactly the same as that for 
Echo II without photoemission. The resulting torque for Echo II with photo- 
emission for orientation I is shown in figure 20(a). 

For orientation II the solar radiation is incident along the Y-axis of the 
satellite coordinate system. Therefore, equation (6k-) must be modified to read 


ai n (P,e) = 


(cos p sln 6 oos ' sln p cos e)ei 
- n a e j sin 0 cos sin a + cos a cos 9 I 


+ 5 x 10"5sin 8 sin t n ep 


A^r^sin 9 d9 


( 77 ) 


The induction torque for Echo II with photoemission for orientation II is given 
in figure 20(b). 

The induction torque for Echo II for the orientation shown in figure 12 
should be calculated next. This induction-torque calculation is a much more 
arduous task than the calculation for the orientation shown in figure 11 inas- 
much as the torques must be averaged at each satellite location over one com- 
plete revolution of the satellite about its spin axis. Because of the more 
difficult calculation and, also, because approximate calculations made by the 
author have shown that there is no appreciable difference in the induction 
torque for the two orientations, the induction torque is not calculated for the 
orientation shown in figure 12. 
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OTHER POSSIBLE SOURCES OF TORQUE 


Several phenomena other than the aforementioned ones may produce torques 
on the satellite. Torques may result from the density gradient in the atmos- 
phere or from the radiation pressure. The resulting torques, however, are 
Believed to he small compared with the eddy-current and the induction torques. 
The aerodynamic torque resulting from the spin should he calculated. Under the 
assumption that all neutral particles striking the satellite are diffusely 
reflected, the aerodynamic torque can he approximated hy an expression given in 
reference 26 as (refer to figs. 4 and 5) 

T z ~ | cos Pi (78) 


For values of p near m/2, equation ( 78 ) is no longer valid and the torque is 
given hy (ref. 26 ) 


T z * |rtpr s 5a3 2 (79) 

The result for these approximations of the aerodynamic torque as applied to 
Echo II is shown in figure 21. In equation (79) the thermal motion of the ions 
has heen neglected. The free molecule drag torque due to the thermal motion is 
given hy Harbour and Lord (ref. 27) as 


T z = ^r s ^nao ^2itmkT = ^ X 10“® newton-meter 


( 80 ) 


The aerodynamic torque is quite small compared with the induction and the eddy- 
current torques. 


ELECTRIC CURRENTS IN SATELLITE SURFACE 


The currents flowing in the surface of Echo II are due primarily to eddy 
currents, photoemission, and accretion of electrons and ions. To compare the 
various currents, a colatitute (3 of O .7 radian was chosen. The eddy current 
(case l) flowing along a gore is, hy equation ( 50 ), 

ig = -j^n:r s 2 B sin a cos T|r n 

This current is constant along a circular ring made up of two gores. For 
P = 0.7 radian the currents in various gores are 
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ig - -b.2 mi 13 i amperes 

for 

0 

11 

d 

ig = -2.4 mi Hi amperes 

for 

f n = 1.00 radian 

ig = - 1.0 milliamperes 

for 

i|r n = l.l4 radians 

ig = - 2.7 milliamperes 

for 

\|r n = 2.43 radians 


The corresponding induction currents without photoemission as determined from 
equation ( 66 ) are shown in figure 22. The induction current is not constant 
along a gore since it is determined hy the accretion of electrons and ions hy 
the satellite. The variation of the induction current for several gores with 
photoemission for orientation I and for orientation II is shown in figure 23- 

The induction current densities for the uniformly conducting satellite 
without photoemission for several values of 0' are shown in figure 24. For 
the uniformly conducting satellite two components of the surface current 
density must he determined since the current is no longer restricted to flow 
along the gores. Also, the induction current densities shown in figures 2b to 
26 are referred to the satellite coordinate system shown in figure l4. 

The variation of induction current density for the uniformly conducting 
satellite with photoemission for orientation I and for orientation II is shown 
in figures 25 and 26 , respectively. 


COMPARISON OF TORQUES 


When comparing the torques it should he noted that all the torques are 
referred to the satellite coordinate system in figure 5 * The results of the 
calculations indicate that only the eddy-current and the induction torques need 
he considered in attempting to find a balance of the torques acting on a satel- 
lite. The induction torque does not depend on the value of the surface resis- 
tivity S used in the calculations. For the uniformly conducting satellite, 

S does not even appear in the equations. For Echo II, D(p) and 0 n (p) are 
hoth proportional to S and their difference is then divided hy R, which is 
also proportional to S. By contrast, the eddy-current torque does depend 
on surface resistivity. However, if some of the connections between the 
connecting ring and the gores represent open circuits, then the value of the 
induction torque will change. 

The average eddy-current torque for the uniformly conducting satellite 
over one earth orbit at co = O .063 radian per second is given in figure 8 as 
^T z ^ = 1.4 X 10 "5 newton-meter. The average torque for arbitrary cd is there- 
fore ^T z ^ = -2.2 X 10”^ru newton-meter. During at least one-half of the orbit 
the satellite will be exposed to sunlight. The value of the induction torque 
will thus be an average of the torques given in figures 17 (a) and 17 (b) - that 

28 


I 



is, (T z ) = -2.2 X 10 - 5 newt on -meter. The induction torque will "balance the 
eddy-current torque for u> = 10 “5 radian per second or for a period of about 
6280 seconds. Thus, for a uniformly conducting satellite of the same size and 
with the same orbit as Echo II, the equilibrium spin rate should be much less 
than that for Echo II. If the satellite is initially not spinning, then its 
rate of spin as determined by the eddy-current and induction torques will be 
given by 



The time constant involved in equation (8l) is about 2 months. 

The average eddy-current torque for Echo II over one earth orbit for the 
orientation shown in figure 10 is found from figure 11 to be 8.6 X 10~5 newton- 
meter. For the orientation shown in figure 12 the average eddy-current torque 
from figure 13 is X 10~5 newton-meter. The actual eddy-current torque for 
Echo II should be bounded by these two values. For the induction torque the 
average of the values given in figures 20(a) and 20(b) - that is 
-3*7 X 10 “5 newton-meter - should be used. It will be noted that the eddy- 
current torque and the induction torque are of the same order of magnitude. 
Thus, considering the uncertainties in the electrical properties of Echo II, 
the agreement between the accelerating and retarding torques is good. 

The exact value of the induction torque will depend on which gores are 
electrically connected to the connection ring and which are not. For example, 
if every other gore is not connected to the top connecting ring (see fig. 9 ) > 
then for a value of (B = 5*97 radians the induction torques are 

Without photoemission: 


T z = -5.8 X 10 “5 newton -meter 
With photoemission, orientation II: 

T z = -1.5 X 10“^ newton-meter 

For the results presented in figures 19 and 20, the corresponding torques under 
conditions where all gores are electrically connected to the connecting ring 
are 

Without photoemission: 


T z = -1.3 X 10~5 newton-meter 


With photoemission, orientation II: 


T z = -1.1 X 10 newton-meter 
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However, when the average torque over one earth orbit is computed for the satel- 
lite with open circuits as previously defined, only a slight increase in the 
induction torque is obtained. It is the author's opinion that any other 
arrangement of open circuits will not appreciably increase the average induc- 
tion torque. Therefore, when electrical continuity is considered, the eddy- 
current torque and the induction torque are found to be nearly equal in magni- 
tude. An interesting factor in the spin history of Echo II is that the spin 
rate increases slightly during the periods when the satellite orbit is 100 per- 
cent in sunlight. These observations are in agreement with the results pre- 
sented in figures 19 and 20 where the induction torque is much larger in the 
presence of photoemission. 

It is of interest to find the time in which Echo II would achieve a steady 
spin rate if <n were zero initially. Note that since the induction torque is 
negative for Echo II, the final also will be negative. For the results 
given in figure 11 (case l), the eddy-current torque for arbitrary <n is 
1.^- X 10“3n newton-meter. Therefore, the spin rate of Echo II as determined by 
the eddy-current and induction torques is given by 

cd = -2.7 X 10“ 2 jl - exp (-1.75 X 10“ 8 t)j (82) 

The time constant involved in equation (82) is about 22 months. 

If the eddy-current torque for the orientation shown in figure 12 
(case II) is used and if the average induction torque is assumed to remain 
unchanged from the previous orientation, then equation ( 82 ) becomes 

cd = -5.3 x 10“ 2 jl - exp (- 8. 8 X 10"9tjj ( 83 ) 

The time constant involved in equation ( 83 ) is about Ml- months. 

Table II summarizes the different torques for Echo II and for the uni- 
formly conducting satellite at the rate of spin indicated. 


CONCLUDING REMARKS 


It was found that only the eddy-current torque and the induction torque 
had to be considered in obtaining a balance of the accelerating and the 
retarding torques. For Echo II ( 196 ^ 4a) the calculated equilibrium rate of 
spin is between -2.7 X 10 “2 and -5«3 X 10-2 radian per second which is close to 
the observed spin rate of - 6.3 X 10 -2 radian per second. Therefore, the prob- 
able explanation of the constancy of the rate of spin of Echo II is that the 
eddy-current torque and the induction torque were in balance and that all other 
torques were negligible. Also, since the induction torque is negative, the 
spin axis should point in the negative z-direction. 
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It should he noted that Echo II was spinning at a rate of about 1 revolu- 
tion per 100 seconds for 9 months after launching. It is also known that this 
spin rate was established during the first orbit. In the light of the present 
analysis, however, it is hard to escape the conclusion that the initial spin 
vector only by chance closely approximated the equilibrium spin vector. 

The surface current was found to be of the order of several milliamperes 
per meter. Whether currents of this magnitude would interfere with experiments 
onboard a satellite will depend on the particular experiment to be performed. 

Langley Research Center, 

Rational Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 1, 1965* 
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APPENDIX 


COMPARISON OF EFFECTIVE SATELLITE RADII 


By use of the method given by Hohl and Wood in reference 10 the potential 
distribution and the ion trajectories have been calculated for two values of |3 
(without photoemission) . The potential distribution is shown in figures 27(a) 
and 27(b) for 0 = 0.4 and 1.39 radians, respectively. The corresponding ion 
trajectories in the plane defined by the satellite velocity vector and the 
Y-axis of the satellite coordinate system (fig. 5) are shown in figures 28(a) 
and 28 (b) . 

Considering the trajectories given in figure 28(a) values obtained for S 
are 7 and 0 Debye lengths for 0=0 and 0 = it, respectively. The corre- 
sponding values obtained by using the result of Jastrow and Pearse (eq. ( 9 )) 
are 5*5 and 0 Debye lengths. For 0 = 1.39* values of 5 obtained from fig- 
ure 28 (b) are 12 and -1.1 Debye lengths, whereas the corresponding values 
obtained from equation ( 9 ) are 15 and -2 Debye lengths. 
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TABLE I.- ATMOSPHERE AND SAT ELLI TE PARAMETERS 


AT 1200-KILOMETER ALTITUDE 


Atmosphere: 

Total mass density, kg/m5 

Ion mass density, kg/m5 

Ion or electron number density, m~5 . . . 

Ion temperature, °K 

Electron temperature, °K 

Debye length, m 

Ion (He + ) thermal velocity, m/sec .... 
Electron thermal velocity, m/ sec .... 

Ion mean free path, m 

Electron mean free path, m 

Ion radius of gyration, m 

Electron radius of gyration, m 

Revolutions per collision for ions . . . 
Revolutions per collision for electrons . 
Magnetic field (over magnetic pole), Wb/m.2 


. . 10-15 
. . 10" 16 
. . 10 10 
. . 10 5 
1.6 x 10-5 


2.7 X 10" 2 

2-5 X 103 

2.7 x 10 5 

2.6 X 10^ 
5.2 X 10^ 
. . 5*6 

5.1 x 10 -2 
1.1 x 10 5 
. . 10 5 

3.8 x 10"5 


Echo II: 

Velocity, m/sec 7-3 X 10^ 

Radius, 20.5 

Spin rate, rad/sec 6-3 X 10“ 2 

Moment of inertia, kg-m 2 7*96 X 10^ 

Surface resistivity, 2 6.7 X 10”5 

Mass, kg 252.9 

Temperature, °K 313 
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TABLE II.- AVERAGE OF TORQUES OVER ONE EARTH ORBIT (t z ) 


Torque 

Echo II 

^cn = - 6.3 X 10 - 2 rad/sec^) 

Uniformly conducting 
satellite 
= -10~5 rad/sec) 

Eddy-current torque, N-m . . . 

1+.3 X 10"5 to 8.6 X 10~5 

2.2 X 10-5 

Surface charge torque, N-m . . 

3.8 x 10- 8 

6 x 10“25 

Aerodynamic torque, N-m . . . 

8.2 X 10 " 8 

1.3 x 10 -9 

Induction torque, N-m .... 

-3.7 x 10 "5 

-2.2 x 10"5 

Coulomb torque, N-m 

-7 x io~ 7 

-7 x 10-7 
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Figure 10.- Satellite orientation for calculation of eddy-current torque. 


Case I. 
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P, radiai 


Figure 11.- Eddy-current 
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Figure 14.- Coordinate system used for calculating Coulomb torque. 
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Figure l6.- Induction torque for uniformly conducting satellite without photoemission. 
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(t>) Orientation II. 

Figure 20.- Induction torque for Echo II with photoemission 
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(a) Orientation I. 
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(t) Orientation II. 

Figure 23.- Variation of induction current for Echo II for several gores with photoemission. 0 = 0.7 radian 
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Figure 25.- Variation of induction current density for uniformly conducting satellite with photoemission 

for orientation I and 3 = 0.7 radian. 
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Figure 26.- Variation of induction current density for uniformly conducting satellite with photoemission 

for orientation II and 0 = 0-7 radian. 
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